Hydrodynamic scaling from the dynamics of relativistic quantum field theory 
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Hydrodynamic behavior is a general feature of interacting systems with many degrees of freedom 
constrained by conservation laws. To date hydrodynamic scaling in relativistic quantum systems has 
been observed in many high energy settings, from cosmic ray detections to accelerators, with large 
particle multiplicity final states. Here we show first evidence for the emergence of hydrodynamic 
scaling in the dynamics of a relativistic quantum field theory. We consider a simple scalar Xcj> 4 model 
in 1+1 dimensions in the Hartree approximation and study the dynamics of two colliding kinks at 
relativistic speeds as well as the decay of a localized high energy density region. The evolution 
of the energy-momentum tensor determines the dynamical local equation of state and allows the 
measurement of the speed of sound. Hydrodynamic scaling emerges at high local energy densities. 
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Hydrodynamics has been used since the work of Lan- 
dau Q|, to describe the properties of high multiplicity 
final states in high energy particle collisions. While Lan- 
dau's motivations dealt with high energy cosmic rays 
there has been since ample evidence from accelerator ex- 
periments that hydrodynamic scaling (the longitudinal 
velocity v x = x/t ), which implies flat rapidity distribu- 
tions, is the correct approximate kinematical constraint 
for the dynamics of high energy particle collisions. 

These findings imply that hydrodynamic scaling must 
emerge from the dynamics of quantum field theories, if 
the latter are to be correct descriptions of collective be- 
havior in particle physics models. While the applicability 
of quantum field theory in these regimes is not in doubt, 
it has not been demonstrated that hydrodynamic scaling, 
which implies that that the energy density isosurfaces are 
surfaces of constant r 2 = (t 2 — x 2 ), is achieved at suffi- 
ciently high center of mass collision energies. 

Present and future experimental prospects for the 
study of hydrodynamic scaling in high energy experi- 
ments are tremendous. The relativistic heavy ion collider 
(RHIC) is presently producing the highest energy, high- 
est multiplicity hadronic final states ever accessible in a 
controlled environment The Large Hadron Collider 
(LHC) will later produce even more spectacular events. 
The detailed understanding of hydrodynamic flows in 
these experiments constitutes the most promising way 
for the determination of the thermodynamic properties 
of nuclear matter at high temperatures || , viz. its equa- 
tion of state, and the nature of the confinement and chiral 
symmetry breaking transition. The connection between 
hydrodynamic scaling and flat rapidity distributions in 
the context of Landau's hydrodynamical model was first 
discussed in [Q. In the context of "boost invariance" , the 
same scaling law was developed in detail by Bjorken Q. 

Direct field theoretical methods, although still in their 



adolescence, offer much promise for the understanding of 
hydrodynamic scaling and the limits of its applicability. 
Moreover they make accessible regimes where particle co- 
herence is important, which escape Boltzmann particle 
methods. In this letter we show, for the first time, how 
hydrodynamic scaling emerges from the dynamics of a 
simple 1+1 dimensional scalar field theory in the Hartree 
approximation. Our results allow us to map the equation 
of state as a function of space and time, and, under well 
known assumptions, determine the speed of sound Cq. 

To exhibit the ubiquity of hydrodynamic scaling we 
study two different situations: one in which a hot region 
is formed in the wake of the collision of two leading parti- 
cles (kinks) at relativistic velocities, and another simpler 
one where we construct a local energy overdensity which 
is allowed to relax under its own self-consistent evolu- 
tion. To be definite we will be concerned with a scalar 
A</> 4 quantum field theory with Lagrangian density 
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The well known Hartree approximation is the simplest 
nontrivial truncation of the coupled equations for the 
field's correlation functions, which assumes that all con- 
nected correlation functions beyond the second are neg- 
ligible ||. This leads to a dynamical equation for the 
mean field tp = (<p) and the connected two-point func- 
tion. We write the quantum field <f> = tp + ip, where tp 
are fluctuations, (tp) = 0. The equations of motion for tp 
and the 2-point function G(x,y) — (ip(x)ip(y)) then are 
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\tp 2 (x) + 3XG(x, x)] tp(x) = 0, 
3X(tp 2 (x) +G(x,x))] G(x,y) = 0. 



(2) 



To solve the equation for the Green's functions we will 
rely on a complete orthogonal mode basis ipk(x) 
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ip(x) 



E [4v* 



x) + a k tp k (x) 



(3) 



where a\.,ak are creation and annihilation operators 
obeying canonical commutation relations. In terms of 
the mode fields ipk, at zero temperature 



(4) 



The effective mass squared of the propagator, x( x , t) 
must be finite which tells us how to choose the bare mass 
n\. In 1+1 dimensions the self-energy has only a loga- 
rithmic divergence, which is eliminated by a simple mass 
renormalization. We choose 



(A = ±m 2 - 3A / = ± 



m 2 -3AG , (5) 



leading to the existence of two homogeneous stable 



phases, corresponding to xo 



m , <fi 



■ and xo = 2m 2 , 



i 2 /X, i.e. a symmetric and broken symmetry 
phase, respectively. The renormalized equations are 

[n+ x (x)-2\<p 2 (x)] <p{x)=0, 

[U + x{x)]^ k {x,t) = V fc , (6) 

X (x) = ±m 2 +3\tp 2 (x)+3\G R (x,x), (7) 

where the renormalized Gr(x,x) = G(x,x) — Gq. The 
challenge posed by Eqs. d(|{?]) in spatially inhomogeneous 
cases is that we need to solve many partial differential 
equations simultaneously. In a spatial lattice of linear 
size L, the computational effort is of order L 2D , per time 
step, where D is the number of space dimensions. Be- 
cause of this demanding scaling we focus on D = 1. 

We consider two classes of initial conditions: 1) Col- 
liding kinks, in the broken phase where, at t=0 we have 

Tfl 

ip(x, t = 0) = -^=^ kink (a; - x )(p kink {-x - x ), (8) 

(9) 



( Pkxak{ x ) = tanh (mx/\/2\ 



with the kinks initially boosted towards each other at 
velocity v, and 2) a Gaussian shape in the unbroken phase 



tp(x, 0) = ipo exp 



x 

2A 



d t (p(x,0) = 0. 



(10) 



In both cases we adopt at t = a Fourier plane-wave 
mode basis, characteristic of the unperturbed vacuum 



Mx,t) = x j —e^+^K 



= V^+X^- (11) 



The orthonormality of the basis is preserved by the evo- 
lution, Eqs. (H-g) 

To study the hydrodynamic behavior we need to spec- 
ify the operator energy momentum tensor T^ v 



T"" = d»<pd v <f> - L. 



(12) 



Its expectation value, in terms of tp and ipk, is 
(Too) = \{v>tf + \(^f + \ E [Itf I" + + Vk, 

k 

(Tu) = \(<Pt? + + \ E M\" + Iffl - Vk, 



x 2 

^ 12A~ A T : 



(13) 



(Toi) = (T 10 > = vtV* + \ E W^t + > 



where all arguments are at x. The subscripts x, t are 
shorthand for spatial and time derivatives, respectively. 
Too and T\\ contain two different types of ultraviolet di- 
vergent contributions. The first arises from the 1-loop 
integral G(x,x). This logarithmic divergence is removed 
by mass renormalization Eq. (0) . The second divergence 
appears in the kinetic and spatial derivative fluctuation 
terms. This divergence is purely quadratic, and is al- 
ready present in the free field theory in the vacuum sec- 
tor. By comparing the mode sum with a covariant di- 
mensional regularization scheme for the free field theory, 
one deduces that the correct subtraction in the mode sum 
scheme is given by 



\ E o^i 2 + m 2 ] -> \ E o^i 2 + i^i 2 - \ k 



I] • (14) 



In practice we discretize the fields ip(x) and the set 
{ifrj} on a spatial lattice with size N and spacing dx and 
use periodic boundary conditions in space. We choose 
dx = 0.125, TV = 1024, and m 2 = 1, A = 1, h = 1. The 
dynamical equations are solved using a symplectic fourth 
order integrator (with a timestep dt = 0.025). With these 
choices, in the finite volume L = Ndx, the momentum 
k takes a finite number of discrete values k n = 2j j^ L , 
with n = { — — l} and continuum k- integrals 
become sums f dk/(2ir) — > L _1 J2 n - The frequency u)k 
now satisfies a lattice form of the dispersion relation, with 



Lul = k 2 +x 2 , k 2 — -— ^(1 - cosdxk n 
dx 



(15) 



These forms also require that the renormalization of T^„ 
be achieved using appropriate lattice choices. In partic- 
ular we adopt \k n \ = ^/k 2 in ([Tj]). 

We are now ready to address the hydrodynamics of our 
field theory. Landau's simplifying assumption, that per- 
vades hydrodynamic simulations of multi-particle flows in 
high energy experiments, is that they behave collectively 
as a perfect fluid, corresponding to 



(Tfj, v ) = (e + p)u^u v - gTp\ d^T^) = 0, 



(16) 



where u M = 7(1,1;), (7 = l/\/l — v 2 ) is the collective 
fluid velocity and e and p are the comoving energy and 
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pressure densities. The latter are the eigenvalues of the 
energy momentum tensor, and can be obtained from the 
invariance of its trace and determinant e — p = TJf, ep = 
Det|T|. The fluid velocity can be obtained from the form 



Eq. (18) implies that energy density isosurfaces lie on 
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(17) 



A perfect fluid is the limiting hydrodynamic behav- 
ior of a collisionless plasma, where transport coefficients, 
such as viscosities are vanishing. In a mean-field approx- 
imation explicit collisions are neglected, so that the form 
(|l6|) is compatible with our dynamical approximation. 




FIG. 1. Contours of equal energy density in space (horizon- 
tal) and time (vertical), near the collision point of two kinks 
(at the origin of the plot), initially boosted towards each other 
at v = 0.8. The collision is symmetric under spatial reflec- 
tion. We show the region after and to the right of the collision 
point. Red denotes the highest energy density, navy blue the 
lowest. The energy isosurfaces show signs of hydrodynamic 
scaling, following approximate hyperboloids, which are dis- 
torted because of the presence of the emerging kinks. 



The attraction of scaling lies in the fact that the re- 
lation x — ±vt allows for significant simplifications of 
the hydrodynamic equations (|l^), which can then be ex- 
pressed in terms of a single variable and thus become 
ordinary differential equations. These can then be solved 
analytically [Q, generating predictions for the spatio- 
temporal behavior of hydrodynamic quantities such as 
energy density e, temperature or entropy. For example, 
for equations of state where dp/de — Cq, where Co is the 
(constant) speed of sound, one easily derives in ID [Q 



e(x,t)/e = (t/x ) 



-(i+cg) . 



(18) 



where eq, xq are integration constants. The equation of 
state p = CqE can be obtained using simple assump- 
tions about the (hadronic) excitation spectrum ||[|. 



hyperboloids t 2 — x A =const, a property that we can eas- 
ily check in our results, see Figs. [I] and || corresponding 
to initial conditions 1) and 2), respectively. Fig. || shows 
the pressure in space-time for the same situation as in 
Fig. |. 




FIG. 2. Pressure isosurfaces in space-time for the situation 
shown in Fig.[l]. Asymptotically far from the kink trajectories 
the equation of state mimics that of a gas with p = CqE. 
Close to the kinks the pressure gradients indicate space-time 
regions where strong energy flows are imminent, such as the 
area around the kink collision point, at the origin of the plot. 




FIG. 3. e isosurfaces in space-time, for the decay of an ini- 
tial Gaussian shape ([h]), with A = 1 and ipo = 5. The energy 
contained in the initial hot region is about two orders of mag- 
nitude larger than that deposited by the kink collision, shown 
in Figs, mid The red region carries away most of the energy 
in the form of a wave-packet traveling close to the speed of 
light. The energy isosurfaces follow exquisite hyperboloids, 
characteristic of hydrodynamic scaling. 
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Besides the shape of energy isosurfaces hydro-dynamic 
scaling also predicts that the single particle distribu- 
tion functions are flat in the particle rapidity variable 
y = I ln[(E + pu)/(E — p\\)], where E and pii refer to the 
energy and momentum in the direction of the collision of 
an outgoing particle. This result is seen experimentally 
in the central rapidity region B , and can be interpreted 
as arising either from the approximate boost invariancc 
of two highly Lorentz contracted colliding nuclei at high 
energies ||, or from the fact that, in the center of mass 
frame, the Lorentz contracted source for particle produc- 
tion has a negligible longitudinal size, when compared 
to the asymptotic particle's (pion) Compton wavelength 
At sufficiently high deposited energies the two ap- 
proaches lead to the same results. 

Scaling solutions do not preserve global energy conser- 
vation. Thus energy isosurfaces must eventually deviate 
from the scaling hyperbola and join neighboring isosur- 
faces, creating characteristic horn like shapes. This be- 
havior, which can be extracted directly from hydrody- 
namic equations, is also observed in the quantum field 
solutions, see Figs. and §. 

Eq. ( |i~8| ) also allows us a measurement of the speed of 
sound Co- Fig. [| shows the decay in time of the energy 
density of the Gaussian mean-field profile already dis- 
cussed in the context of Fig. ||. The resulting fit suggests 
a value of 1 > cq Z 0.77, compatible with an ultrarela- 
tivistic equation of state cq = 1 in ID. 




FIG. 4. e(x = 0,t), in the wake of the decay of a Gaus- 
sian hot region of Fig. ^. The solid lines show power laws of 
the form ( |l8| ) with 1 + Cq = 2 (blue), corresponding to the 
ultrarelativistic ideal gas (co = 1) and with 1 + c\ — 1.6. 



in a variety of settings, that hydrodynamic scaling 
emerges from the dynamics of quantum field theory at 
sufficiently high energy densities. We analyzed situations 
both with leading particles, which constitute the asymp- 
totic states both before and after the collision, and fol- 
lowing the evolution of a simple local energy overdensity. 
The extension of this type of calculation to 3D, where 
hydrodynamics is richer, and to include scattering, nec- 
essary for the description of real fluids, remain necessary 
steps to make real time studies of quantum fields predic- 
tive experimentally in the context of heavy ion collisions. 
Nevertheless studying the real time dynamics of quantum 
field theories demonstrates the applicability of models of 
particle physics in the largely unexplored limits of very 
large time and spatial scales, where they acquire fasci- 
nating macroscopic and non-perturbative properties. 
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